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In the non-relativistic theory of gravity recently proposed by Hofava, the Hamiltonian constraint 
is not satisfied locally at each point in space. The absence of the local Hamiltonian constraint 
allows the system to have an extra dark-matter-like component as an integration constant. We 
discuss consequences of this fact in the context of cosmological perturbations, paying a particular 
attention to the large scale evolution of the curvature perturbation. The curvature perturbation is 
defined in a gauge invariant manner with this "dark matter" taken into account. We then clarify 
the conditions under which the curvature perturbation is conserved on large scales. This is done by 
using the evolution equations. 

PACS numbers: 04.60.-m, 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

A power-counting renormalizable theory of gravity, 
proposed recently by Hofava 0, 0] , has attracted much 
attention. The essential aspect of the theory is broken 
Lorentz invariance in the ultraviolet (UV), where it ex- 
hibits a Lifshitz-like anisotropic scaling, t — > £ z t, x — > Ix, 
with the dynamical critical exponent z = 3. This will 
bring an interesting change in the physics of the early 
universe since the UV effect may play an important role 
there. The study of the cosmology based on Hofava grav- 
ity, which is called Hofava-Lifshitz (HL) cosmology, has 
been initiated by Refs. [HQ, and since then various as- 
pects of HL cosmology have been explored, including the 
generation of chiral gravitational waves [f| , a new mech- 
anism to generate a scale invariant primordial spectrum 
without inflation Q, the bouncing scenario Q, and oth- 
ers Aside from cosmology, other interesting works 
can also be found in 

Several versions of Hofava gravity have been known, 
which arc classified according to whether or not the de- 
tailed balance and the projectability conditions are im- 
posed. Among them the theory with projectability and 
without detailed balance is argued to evade the prob- 
le ms fiol [ll| pointed out in the literature [H, EH, [HI, 
Il5llld| (see also [13] )• The most distinguished feature of 
projectable Hofava gravity is that the Hamiltonian con- 
straint is not a local equation satisfied at each point in 
space, but rather a global equation integrated over the 
whole space. Since the global Hamiltonian constraint is 
less restrictive than the local one, it allows for a wider 
class of solutions which contain an additional dust-like 
component as an integration constant, as was clearly re- 
marked in [Tol | ■ 

In this paper, we discuss consequences of the absence 
of the local Hamiltonian constraint within the context 
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of cosmological perturbations. Cosmological perturba- 
tions in Hofava gravity have already been investigated 
in [H, 03 HI El, [H, HI, IH, HI, [H, many of which 
are interested in the modified dispersion relation in the 
UV regime. Our focus here is on the large scale evolu- 
tion of the curvature perturbation in the Sotiriou-Visser- 
Weinfurtner generalization of Hofava gravity [13] ■ Wang 
and Maartens were the first to study cosmological per- 
turbations in this version of the theory (25[. In general 
relativity, the curvature perturbation on uniform den- 
sity hypersurfaces, commonly denoted as Q, is conserved 
on large scales, provided that the non-adiabatic pressure 
perturbation is negligible. This fact can be proven by 
utilizing the energy conservation law only [28| , and hence 
the conservation of £ holds true in a wide range of grav- 
ity theories such as brane- world gravity [29[. In Hofava 
gravity, however, an extra degree of freedom mimicking 
dark matter, or, what is dubbed as "dark matter as an 
integration constant" in [l(| [H|, appears as a natural 
outcome of the lack of the local Hamiltonian constraint. 
This forces one to consider effectively a multi-fluid sys- 
tem even if the system is composed of a single (real) fluid, 
which implies that the effect of the entropy perturbation 
may not be negligible. Furthermore, individual compo- 
nents are in general not conserved separately, while the 
total energy including "dark matter" thus introduced is 
shown to be conserved locally (by invoking the evolution 
equations). Therefore, we start with defining the gauge 
invariant curvature perturbation, emphasizing the pro- 
jectability condition and taking into account the presence 
of "dark matter as an integration constant." Then, we 
discuss the conditions under which the curvature pertur- 
bation is conserved on large scales by using the evolution 
equations. 

This paper is organized as follows. In the next section, 
the basic equations in Hofava's non-relativistic theory of 
gravity is provided. We review the background cosmol- 
ogy in Hofava gravity in Sec. Ill, emphasizing the conse- 
quence of the absence of the local Hamiltonian constraint. 
Then, in Sec. IV, we study the large scale evolution of the 
curvature perturbation in the absence of the local Hamil- 
tonian constraint and clarify the conditions under which 
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it is conserved on large scales. We draw our conclusions 
in Sec. V. 



II. HORAVA GRAVITY 

We consider the projectable version of Hofava gravity 
without detailed balance. The dynamical variables are 
M, Mi, and 7y , in terms of which the metric can be 
written as 



d S ^ 



M 2 dt 2 + 7 y- (dx i + M l dt) (dx 3 + M j dt) , (1) 



with M l = j^JVj. The projectability condition states 
that the lapse function depends only on the time coordi- 
nate, M = M(t), while Mi and jij may depend on t and 
x. The theory is invariant under the foliation-preserving 
diffeomorphism: t — > t(t), x l — > x l (t,x). Under the in- 
finitesimal transformation. 



t^t + X °(t), x^x*+ X \t,x), 
these variables transform as 



(2) 



lij Hj-jijX -lik^jX -IjkViX , 

M -> M-M x °-Mx a , (3) 
Mi ^ Mi- W iX j Mj - x^jMi - x j Hj - x°Mi - X °M, 

where a dot is the derivative with respect to the time 
coordinate t and Vj the covariant derivative associated 
with the spatial metric %j . One can see that M remains 
x-independent after the transformation, and thus it is 
natural to impose the projectability condition. 

The dynamical variables are subject to the action 0, 

S3 

S = J dtd 3 xM^f (KijKV - XK 2 + R + C V2 ) 

+ J dtd^xM^Cm, (4) 



where C m is the Lagrangian for matter fields, 



(5) 



is the extrinsic curvature, R = j^Rij is the trace of the 
Ricci scalar (the spatial curvature scalar), and 

C V2 := a 2 R 2 + a 3 R i3 R i3 + a 4 R 3 + a b RR i3 R i3 
• .»,./.',' /.'.'/.V + a 7 RV t V l R + osV./fj/V //'*'((>) 

The kinetic term coincides with that of general relativ- 
ity when A = 1, but we do not specify the value of A 
throughout the paper. Note however that it has been 
argued that an additional longitudinal degree of free- 
dom of gravitons suffers from ghost-like instabilities for 
1/3 < A < 1 @ , . One may include a cosmological con- 
stant in the above action, but we do not write it explicitly 



since it can also be thought of as a part of the matter La- 
grangian. One can also include a parity-violating term 
associated with the Cotton tensor, as in original Hofava 
gravity 0). 

Variation with respect to M yields the Hamiltonian 
constraint. In the projectable version of Hofava gravity, 
the Hamiltonian constraint is not satisfied locally at each 
spatial point, but rather a global equation integrated over 
the whole space because M is a function of t only. The 
global Hamiltonian constraint reads 

J d 3 x^ [KijK ij - \K 2 -R- Cv-2 + l^GE] = 0,(7) 



where 



-C m - M 



SM 



(8) 



Variation with respect to Mi leads to the momentum con- 
straint, 



where 



V,/"' = 8ttGJ\ 



pij .- K ij _ XKj ij , J* = -M 



SMt 



(9) 



(10) 



Finally, variation with respect to 7 jj gives the evolution 
equations, 

2 {K ih K* - XKKi 3 ) - i [K u K kl - XK 2 ) 7y - 



My/j 



-lv fc (PikMj) + -iv fc (P jk Mi) + Ra - \Ri 



+Fij = 8irGT ij} 
where := 5£v2/S r y lj - (l/2)~/ijCv2 and 

. 5 C m 



Tij . £ m ^fij 2 



2" i*3 

^ (11) 



[12) 



The matter action is invariant under the infinites- 
imal transformation (J3J), which results in the energy- 
momentum conservation equations: 



d 6 x 



0(13) 



~{V j M i -V i M j ) = Q. (14) 

The energy conservation law (|13p is of the form of the 
integration over the whole space, as is the case for the 
Hamiltonian constraint. 
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III. BACKGROUND EVOLUTION 

The background evolution of HL cosmology can be de- 
rived by setting M = 1, Mi = 0, and 7^ = a 2 (t)Sij. The 
evolution equation at zeroth order reads 



1 - 3A 



3H 2 



2H) = 8nGp, 



(15) 



where H := a/a, and Xy = pr/ij has been assumed. 
Let us define £(t) by 



8nG[£{t) + p\ 



°-{l^3X)H 2 , 



(16) 



where p is the background value of the matter energy 
density E. In the case of A = 1, the meaning of £ becomes 
more transparent by noticing that 8irG£ = 8ttGT ° — G ° , 
where G ° is the (0 0) component of the usual Einstein 
tensor: £ arises because the local Hamiltonian constraint 
is absent in Hofava gravity. This term corresponds to 
"dark matter as an integration constant" in Refs. (Iol[Tl|. 

We emphasize that in this paper the homogeneous 
background is assumed at least in our observable patch 
of the universe because nobody can tell what happens 
beyond the present horizon scale. Under this assumption 
we can conclude from the global Hamiltonian constraint 
that £ does not necessarily vanish in the local patch. For 
example, we may have £ > in our patch of the universe, 
but £ may be negative in a different patch. Our assump- 
tion is in contrast to Ref. [25| , in which the standard as- 
sumption of a homogeneous background is made over the 
whole space, so that the global Hamiltonian constraint 
enforces £ = 0. In this paper, we do not assume the ho- 
mogeneity over the whole space, and therefore, Eq. ([7]) 
does not constrain the value of £ in our observable patch 
of the universe. 

In terms of £ , Eq. (|15p can be written in the form of a 
conservation equation: 



£ + p + 3H(£ + p + p)=0. 



(17) 



This does not guarantee the local conservation of the 
matter energy density. If the matter action respects gen- 
eral covariance, we have an additional conservation equa- 
tion, p + 3H(p + p) = 0. In the case of scalar field mat- 
ter [J], the equation of motion leads to p + 3H(p+p) = 0. 
Combining the local conservation of matter energy with 
Eq. (fTT|) , we obtain £ + 3H£ = (^implying that £ indeed 
shows a dust-like behavior [h], [ll[ . 



IV. LARGE SCALE COSMOLOGICAL 
PERTURBATIONS 

Let us study linear perturbations around the cosmo- 
logical background. The perturbed metric is given by 



M 2 = 1 + 2A(t), Mi = a 2 B^ 
7ij = a 2 [(1 - 24,)S. lJ + -21) , , 



(18) 



Since we are imposing the projcctability condition, the 
perturbation of the lapse function A does not depend 
on x. Cosmological perturbation theory in Hofava grav- 
ity without the projectability condition has been studied 
in [2l[. Under the scalar gauge transformation, i.e., the 
infinitesimal transformation with x % = d l x(t,x), the 
metric perturbations transform as 



A^A-x°, yj^ip + H x °, 
B —> B — x, D —> D — x- 



(19) 



Since x° depends only on t, inhomogencous ip cannot be 
gauged away, while A can be set to zero by the gauge 
transformation. This point is in contrast to general rel- 
ativity. It is convenient to define a := D — B, which is 
gauge invariant. 

The evolution equations take the form 



gs/ + ^ - 

= 8irG6T/, 



1 



& + 3Ha 



SF/ 
(20) 



where 



G := -(1-3A) 



i> + 3H*f} + HA+ 3H 2 + 2H) A 



-(1 - A)V 2 [& + 3Hcr] - -V s 



& + 3Ha 



(21) 



and SFj 3 is to be derived from Cv2- Since SF^ = 
©(V 4 ), 1 it is not important as long as one concerns the 
large scale evolution of cosmological perturbations. (We 
do not consider the case in which higher spatial derivative 
terms are much larger than C(V 2 ) terms, though C(V 6 ) 
terms have an interesting effect on the spectrum of per- 
turbations @.) The perturbed energy- momentum tensor 
may be written in terms of isotropic and anisotropic pres- 
sure perturbations as 



5p 8? + di&> 



1 



-V z 8, 



2a i 



n. 



(22) 



The momentum constraint is given by 



d i -(1-3A)?A + (1-A)V 2 ( t =8nGa 2 SJ l . (23) 

Analogously to £ defined in the previous section, let 
us now define e(t, x) by 



8irG[e{t,x) +Sp] = -3(1 - 3X)H (ip + HA 



f 4 + Ha 



(24) 



A straightforward calculation shows that Rij = didjip + V 2 i/>i5ij 
and the variable D does not appear here. Therefore, D 
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where Sp is the perturbation of the matter energy density 
E given in Eq. ©. In the case of A = 1 we have 8wGe = 
8nGST ° — SGq, from which it is clear again that e is 
a consequence of the absence of the local Hamiltonian 
constraint. Thus, e may be regarded as a energy density 
perturbation of "dark matter as an integration constant." 
It is easy to check that e transforms as e — > e — £x° under 
the gauge transformation. In terms of e, the trace part 
of the evolution equations, Q — 8nGSp, can be written 
as 

i + 5p + 3H (e + 5p + Sp) - 3ip (£ + p + p) 



1 



8-kG 
+o (V 4 ' 



1> 



2Ha + 3H(l- A) {a + 3Ha) 



(25) 



Sp + 3H (Sp + Sp) - 3^{p + p) = 0(V 2 ). 



which reminds us of the perturbed energy conservation 
equation. If the matter energy density is conserved locally 
at perturbative order on large scales, one has, in addition 
to Eq. HSJ, 



(26) 



Note, however, that Eq. (|25|) does not necessarily imply 
the local conservation of the matter energy density (|2"6")) . 

In the following we shall study two different cases £ 
and £ = 0. 



A. £ ^0 

Formally, one can define the following gauge invariant 
quantities: 



C := (1 - /)Chl + /C„ 



(27) 



Chl:=-V-^, Cm-=-1>-H-Z, (28) 
£ P 



with 



/(*) 



£ + P 



(29) 



Only when the matter energy density is conserved locally 
at zeroth order, we may rewrite Eq. (|29p to have 



/ 



p + p 



(30) 



£ + p + p 

Using the definition (f2"4"|) , ( can be written more explicitly 



as 



H \ > 3 1 



2V 2 (tp/a 2 + Ha) 
{1-3X)H 



(31) 



so that C can be expressed solely in terms of the metric 
perturbations. This is essentially the same as the quan- 
tity first introduced in [3(|. Using the variables defined 



above and neglecting the C(V 2 ) terms, we can rewrite 
Eq. ([23]) in a suggestive form as 



c 



H 



P 



-Sp 



nad 



Hc 2 J(l - /)5 HL , (32) 



where we have introduced the non-adiabatic pressure per- 
turbation of matter, <5p n ad := 5p — c 2 s 8p, with c 2 := p/p, 
and the isocurvature fluctuation between "dark matter 
as an integration constant" and ordinary matter, 



Shl '■= 3 (Chl — Cm) 



(33) 



We emphasize that Eq. (|32j) has been derived only by 
using the evolution equations. 

Equation (|3"2")) however tells nothing about the large 
scale evolution of C unless the evolution of Cm and Chl 
is specified (except for the special case c 2 /(l — /) ~ 0). 
If the matter energy is conserved locally, one finds, from 
Eqs. (QUI and (J2BJ), that / + 3Hc 2 J(l - f) = at zeroth 
order and Cm — on large scales at perturbative order, 
assuming that <5p n ad = 0. In this case it can be shown 
that Chl — on large scales, but C is n °t conserved in 
general. Indeed, it follows immediately from the defini- 
tion Eq. (@ZD that ((t,x) = [1 - f(t)](£l(x) + f(t)(£ ] (x), 

where Chl an d Cm^ are the initial conditions for the cor- 
responding variables. For dust-like matter with p oc a~ 3 , 
/ is constant since £ also scales as a~ 3 , and hence C is 
conserved. This fact was already clear in Eq. ((32)) with 
c 2 = 0. Another case in which C is conserved on large 
scales is 6>hl = 0, that is, Cm = Chl- Whether this "adi- 
abatic relation" between "dark matter" and usual mat- 
ter is likely or not depends upon the specific scenario in 
the early universe. For example, in the case discussed 
in Ref. [y], ordinary matter is produced by the decay of 
the curvaton or the modulus while the initial condition 
of "dark matter" is determined by quantum fluctuations 
of the combination of scalar gravitons and (real) matter 
fields. In this case, there is no relation between them in 
general, leading to Shl 0- 

Interesting cases with / ~ and / ~ 1 can be studied 
without knowing the evolution of <Shl and hence without 
relying on the local conservation of the ordinary matter 
energy density. If radiation dominates the energy den- 
sity of the universe at early times, we have (1 — /) ~ 0, 
which leads to the conservation of C during that period. 
On the other hand, if "dark matter as an integration 
constant" accounts for a significant portion of real dark 
matter and dominates the energy density of the uni- 
verse at late times, we have / ~ 0, which again leads 
to the conservation of C- However, in the intermediate 
regime, /(l — /) = 0(1), so that the curvature pertur- 
bation grows provided that Skl ^ 0. This property is in 
accordance with what is found in a conventional multi- 
fluid system [28| . One should also notice that in the case 
where "dark matter as an integration constant" consti- 
tutes a large portion of real dark matter, 5hl represents 
the isocurvature fluctuation between dark matter and ra- 
diation, which is strongly constrained by the cosmic mi- 
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crowave background (CMB) anisotropy. Thus, the sce- 
nario in which "dark matter as an integration constant" 
is really a dark matter component and there is no natural 
reason to explain 6>hl = gives rise to a large isocurva- 
ture fluctuation, which could be incompatible with the 
present constraint. 

Since A does not depend on x, taking the spatial gra- 
dient of Eq. (jni) yields 



di 



H 



H 



(34) 



(One can do essentially the same thing by making the 
gauge choice ^4 = instead.) This equation is useful for 
reconstructing the curvature perturbation from £. 
Note that ditp is gauge invariant because tp is subject only 
to the temporal gauge transformation t — > t+x°(t). lidi( 
is constant in time, one finds diip = —di((x) + H(t)Ci(x), 
where the second term corresponds to the decaying mode. 

For an illustrative purpose let us consider conserved 
matter with the equation of state p = wp. Equation 
is integrated to give 



diip(t, x) 



3(1 + w)„,. s f* At 1 



-H{t) 



i + s/p- 



(35) 



For £ <C p, we obtain diip ' 

Substituting the curvature perturbation diip into the 
traceless part of the evolution equations ([20)) , we obtain 
the metric shear di<r. Once these metric perturbations 
are determined, the (large scale) CMB anisotropics can 
be computed by making use of the perturbed geodesic 
equations. The detailed calculation of the CMB tem- 
perature anisotropics in HL cosmology is left for further 
study. 



B. £ = 

In this case it does not make sense to define £hl, but 
still one may define ( directly by Eq. (|3"Tj) . The trace 
part of the evolution equations implies 



H 



P 



(Spna.d 



(36) 



(Note that s is gauge invariant when £ = 0.) Thus, in 
general £ is not conserved even if Sp na d = 0. 

Let us consider again the simple case where the matter 
energy is conserved locally and 5p na( \ = 0. In this case it 
is easy to see that e = Eo/a 3 . Integrating Eq. (|3"6"1) . we 
obtain 



(0) 



Kp + p) 

where to is some initial time 



Hp + pY 



(37) 



V. CONCLUSIONS 

In this paper, we have studied the large scale evolution 
of the cosmological curvature perturbation in projectable 
Hof ava gravity, emphasizing the effect of "dark matter as 
an integration constant" [ill 03 that appears as a conse- 
quence of the global Hamiltonian constraint. Our view is 
that we cannot tell the cosmological dynamics far outside 
the present Hubble horizon, and hence the global Hamil- 
tonian constraint does not provide any information in 
our observable patch of the universe. This assumption 
makes the impact of "dark matter as an integration con- 
stant" rather non-trivial. The curvature perturbation £ 
has been defined in a gauge invariant manner with this 
"dark matter" component taken into account. We then 
clarified the conditions under which £ is conserved on 
large scales by invoking the evolution equations. In par- 
ticular, we pointed out that £ is sourced by the relative 
entropy perturbation 5hl between "dark matter as an 
integration constant" and ordinary matter. This source 
term is effective during the period when c 2 s f{l — /) is not 
negligible. In that period, we need to know the evolution 
of 5hl in order to know the evolution of C- This is made 
possible by assuming the local conservation of the energy 
density of ordinary matter. 

If the "dark matter" component constitutes a large 
portion of real dark matter, Shl corresponds to the 
isocurvaturc fluctuation between radiation and dark mat- 
ter, which is strongly constrained by the cosmic mi- 
crowave background anisotropy. In this case, one there- 
fore needs a natural reason to explain <Shl — 0, which is 
a challenge in HL cosmology. 

In the present paper, we have focused on the supcrhori- 
zon evolution of the curvature perturbation for a given 
initial condition. In order to impose an appropriate ini- 
tial condition, we need to specify the scenario of the early 
stage of the universe and then to quantize the cosmolog- 
ical perturbations, which would allow us to give obser- 
vational prediction for curvature and isocurvature per- 
turbations. Although the procedure is familiar and es- 
tablished in conventional cosmology, it is non-trivial in 
HL cosmology. It would be interesting to study quanti- 
zation of the coupled system of a scalar field and metric 
perturbations and solve its evolution from subhorizon (or 
WKB) to superhorizon regimes in HL cosmology. More- 
over, in HL cosmology we have a novel mechanism to 
generate a scale invariant spectrum of quantum fluctua- 
tions [|| , which relies on the modified dispersion relation 
brought by C(V 6 ) terms and does not require inflation. 
A detailed analysis of this mechanism taking into account 
the effect of the metric perturbations is yet to be done. 
These issues are left for a future study. 
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